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1.  Introduction 


Mountain-valley  wind  circulation  has  long  been  of  Interest  to 
Inhabitants  in  mountainous  areas  and  meteorologists.  Yet  it  has  not  been 
fully  explored  and  xuiderstood.  Various  daily  operations  require  much  more 
detailed  meteorological  knowledge  than  in  the  past.  In  particular,  an 
accurate  short-range  small  meso-scale  3-diraensional  prediction  model  is 
needed  over  mountain- valley  terrain.  Aiming  at  this  goal  we  shall  divide 
the  investigation  into  three  tasks;  l)tp  develop  a 3-dimensional  numerical 
model  with  the  initial  value  approach  and  to  study  the  basic  physics  in 
this  particular  kind  of  local  wind  system,  2)  to  compare  results  from  the 
numerical  model  with  some  simple  analytical  solutions  and  available  obser- 
vations, and  3)  to  build  a short  range  meso-scale  prediction  model  for 
mountain- valley  wind  systems. 

This  report  discusses  problems  associated  with  the  first  task  and 
partially  with  the  second.  The  last  task  and  part  Oa  the  second  one  shall 
be  studied  as  a continuation  of  the  tasks  presently  undertaken  and  will  be 
completed  in  the  next  phase  of  investigation. 

Theoretical  and  numerical  studies  of  flow  over  a small  mountain  ridge 
have  been  numerous.  The  most  recent  studies  by  Mahrer  and  Pielke  (1975.1977) 
and  Clark  are  interesting.  However,  the  boundary  layer  problem,  such  as 
flow  in  a typical  valley  region,  needs  further  attention  and  investigation. 
Both  theoretical  and  numerical  sutdies  of  mountain-valley  wind  systems  for 
a typical  3-diMQ‘Tii'ional  meso-scale  mountain-valley  terrain  have  bean  seldom 
found  in  the  literature,  i-ew  intrigue  linear  or  non-linear  interactions 
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between  the  free  atmospheric  flow  and  the  thermally  induced  flow  near  the 
surface  in  a typical  valley  have  been  reported  except  for  the  work  by 
Tang(1976),  Tang  and  Peng(l97^)  for  2-dimensional  ones,  and  over  general 
mountain-basin  situation  by  Kahrer  and  Pielke(1977)t  Studies  on  the  inter- 
play between  the  slope  wind  or  the  cross-valley  flow  and  the  valley  wind 
(along  the  valley  axis)  are  completely  lacking  in  the  literature. 

In  studying  mountain-valley  wind  systems  the  non-hydrostatic  approxi- 
mation is  believed  to  be  important  because  of  the  basic  convective  nature 
of  the  problem,  especially  when  horizontal  and  vertical  scales  become  com- 
parable. Since  this  investigation  emphasizes  the  study  of  the  flow  in  the 
boundary  layer,  computations  with  parameterizations  under  various  meteoro- 
logical conditions  for  the  present  model  have  been  performed.  The  develop- 
ment of  the  numerical  model  is  presenter  in  Section  II.  Resvilts  of  the 
numerical  modeling  under  various  conditions  are  presented  in  Section  III. 
Some  available  observational  data  will  be  used  for  comparison  and  discussed 
at  the  end  of  Section  III.  Discussion  of  the  study  and  future  study  will 
be  in  Section  IV. 


II.  The  development  oi'  numerical  model 
1)  The  differential  equations 

In  order  to  simulate  circulation  over  a valley  terrain,  the  configura- 
tion of  the  relief  and  the  slope  of  the  topography  are  believed  to  be  critical. 
In  contrast  to  the  study  of  lee  waves  for  which  linearized  lower  boundary 
surface  condition  may  be  used  to  obtain  a general  picture  of  the  wave  pattern 
in  the  middle  of  the  atmosphere,  non-linear  lower  boundary  conditon  should  be 
used  here.  Furthermore,  a set  of  equations  should  have  coordinates  trans- 
formed so  that  it  can  be  used  for  arbitrary  terrain.  A generalized  <r-  . 
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coordinates  In  z Is  taus  used  In  the  model.  The  basic  transformation  is 


H TT-S-  (1) 

where  z = height  in  rectangular  coordinate  system, 
z's  height  in  transformed  coordinate  system, 
h(x,y)=terrain  height  profile  in  regular  rectangular  (x,y,z)  coordinates, 
and  H = height  of  top  boundary  over  horizontal  ground  for  the  model. 


With  this  transformation  we  obtain  the  following  relationships; 


1 

\ H-f) 


^ St 


Note  that  the  square  brackets  with  subscripts  have  special  meanings  in  the  equations, 

Similar  results  can  be  obtained  for  the  differential  with  respect  to  y.  A 
set  of  governing  differential  equations  for  an  incompressible  atmosphere 
with  Boussinesq  approximation,  after  chopping  off  the  prime  for  transformed 
cooz*dinates,  can  be  written  as: 

(u)^  ^ _ [Fj^+  -fv  n-4-,«a/  (2) 


3 


r 


[I 


(v)^t  [ kv]  + [v 

■h  v»-(k^  ''V-) 


M 


(3) 


t(  [ h,yi^-4  w (>^K  ( « Y«)^[['«4^'‘]^  ‘^“'•(1^ 

C >^3^-*- Mj.-+ =o 

where  ( )ti  ( )a,  etc  = partial  differential  with  respect  to  t.  and  z.  etc 

respectively  in  transformed  coordinates. 
u,v,w  = velocities  in  the  x-,  y-,  and  z-  axes  directions. 
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and  = pressure  deviation  from  the  mean  at  a given  height 

= ^ ir' ^ TT^if  ) 

2)  The  topography 

The  2-dimen3ional  valley  cross-section  and  mesh  in  the  x-z  plane  at 
y=0  is  shown  in  figure  1.  The  x-axis  oriented  in  the  cross-valley  direction, 
the  y-axis,  in  the  down- valley  direction,  and  the  z-axis,  in  the  vertical 
direction  pointed  upward  from  the  gound  surface.  The  grid  interval  in  the 
vertical  direction  from  ground  up  are  50.50,100,300,500,1000  m up  to  the 
sixth  grid  and  lOOOra  thereafter  for  every  grid  interval  up  to  10km  level. 

The  grid  size  in  the  x-  or  y-axis  direction  is  even  but  is  different  for 
different  computer  runs.  The  grid  size  in  the  x-axis  direction  ranges  from 
333®  to  8km,  and  in  the  y-axis  direction,  from  5km  to  15km  depending  upon 
the  run.  Experiments  with  combinations  of  different  size  grids  in  the  x- 
and  y-axas  directions  are  made  for  purposes  of  comparison.  The  time  interval 
chosen  depends  upon  the  grid  size  used  in  order  to  fulfil  the  Courant- Friedrichs 
-Levy  linear  computational  stability  criterion. 

The  height  of  ridges  above  the  valley  floor  is  assumed  to  bo  constant 
in  the  y-axis  direction,  but  it  varies  from  case  to  case,  ranging  from  100m 
to  1.5km.  The  height  of  valley  floor  at  the  valley  head  is  300m  above  that 
at  the  valley  mouth  at  the  bottom  of  the  valley.  The  cross-section  along 
the  valley  floor  is  shown  in  Fig.  2. 

3)  Boundary  layer 

Since  mountain- valley  systems  are  largely  affected  by  differential 
heating  and  terrain  features,  some  fine  turbulent  characteristics  over 
the  uneven  lower  boundary  are  indeed  import?int  to  the  investigation.  It 
vrould  be  interesting  to  have  some  good  knowledge  of  vertical  and  horizontal 
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eddy  viscosity  and  diffusivity  particularly  suited  to  the  valley  terrain. 
Unfortunately,  little  literature  about  this  in  a valley  is  known.  Because 
of  lack  of  the  precise  knowledge  in  this  matter,  the  eddy  viscosity  and 
diffusivity  over  homogeneous  plane  surface  formulation  may  be  adapted  or 
a simple  constant  with  a relatively  large  numerical  value  is  to  be  used. 

The  relatively  large  value  is  used  because  of  the  basically  strong  turbulent 
nature  over  rough  mountainous  surface.  In  this  report  we  have  undertaken 
the  study  with  both  approaches. 

Adapting  the  basic  formulation  due  to  Businger(  1973) . ''‘0  can  compute 
eddy  viscosity  and  diffusivity  by  solving  well-known  equations.  Iva  first 
determined  the  Richardson  number, R^,  from  the  mean  vdna  speed  at  one  height 
(say  the  first  grid  point  above  ground  surface)  and  temperature  at  two 
heights  (say,  first  grid  level  and  surface).  By  using  the  profile  Richard- 
son number,  ^ ^ may  be  expressed  as 


(i- 

i 


A 

H u.'jiy 


4ori^ 


o 


(7) 


> D 


where  z is  non-dimensional  height  which  is  the  ratio  of  the  height  above 
the  ground,  z,  and  the  Konin-Obukov  length.  Determining  z from  the.  last 
equation  and  therefore  Monin-Obukov  length  from  formula  L = z/  z lead 
to  the  determination  of  non-dimensional  momentum  and  temperature,  ^ and 
(p  respectively,  from  the  following  relationships  with  z, 


^ r I 

1 I -J  ^ 


2^0 


^ > O 
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IL 
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^ o.7<^ 


O.74  0-^zy'''^ 


2 ^ tP 
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where 
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Aa^  2f^ 


Hi 


* ^u  S 

TT'  > 
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Tf 
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The  integrated  version  of  the  profiles  are,  respectively, 

e.itfWi-i'xS 


,t.  = Ail — 

Mp-fi  ' 


9,  = 


(10) 


where 


^ 2 •+  ^ 

{ 

A 

2;  > 0 

(11) 

■-  2 

2,  > CP 

and 


““  A'=fO-'  . -=(<3»A' 

The  values  of  (j)^  and  will  allow  first  determination  of 

through  the  above  set  of  equations  and  then  14.^  and  6>if  with 
known  Zq  and  ^ (zq)*  Oi^ce  and  "tCy  are  determined,  we  determine 

vertical  eddy  viscosity  and  diffusivity,  and  respectively 

in  the  surface  layer  by 


Y-»| 


(12) 


and 


and  trTT-  Ky  by 
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The  determination  of  these  quantities  is  consistent  with  additional  assump- 
tions. Above  the  surface  layer  we  use  O'Brien  's  (1970)  interpolation, 
i.e.  let  K stand  for  or 


= the  height  at  the  top  of  the  planetary  boundary  layer. 
Ky(zi)  will  be  sot  1 cm^seo"''  and  will  be  forecasted  by  an  equation 
proposed  by  Deardorff( 1974)  and  tested  by  Pielke  and  Mahrer( 1975) . 


di' 

where 


J 3 


•I/S'm 


-f 


= the  stability  immediately  above 
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V>/y  = the  vertical  velocity  at 
The  depth  of  the  surface  layer,  will  be  taken  as  0.04  z^. 


(15) 
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4)  The  boundary  and  Initial  conditions 


At  the  top  boundary  where  z = H = 10km,  the  horizontal  velocity  is 
speciiled  as  u = U and  v = 0,  Since  the  top  boundary  is  taken  as  a flat 
solid  surface,  it  is  natural  to  assume 
•3*  3 ^ 

•w  « a,  s —tc"  ss.  o 

At  the  ground  it  is  assumed  that  no  slip  condition  is  used  for  wind 
(i.e,  u=v=w=0).  The  air  temperature  at  surface  can  be  calculated  from  the 
energy  balance  equation  at  the  ground.  In  order  to  apply  the  present  model, 
an  analytical  expression  of  the  ground  surface  temperature  and  heat  flux 
at  each  small  time  step  in  terms  of  solar  and  long  wave  radiation,  geome- 
trical factors  of  the  terrain,  and  initial  surface  temperature  has  been 
derived,  i^etails  of  the  investigation  is  shown  in  Appendix  iA.  In  this 
report,  however,  only  specified  temperature  variatiors  on  the  ground  are 
used.  The  computed  temperature  will  be  incorporated  into  the  present  model 
in  the  next  phase.  The  boundary  conditions  for  J can  be  obtained  from 
the  third  equation  of  motion.  At  lateral  boundaries  all  are  assumed  to  be 
periodic. 

The  initial  condition  is  chosen  such  that  all  quantities  are  zero  in 
space  except  for  the  specified  vertical  distribution  of  moan  U with  height. 

In  order  to  avoid  the  creation  of  relatively  strong  disturbances  at  the 
initial  stage,  we  let  the  mountain  be  flat  in  the  beginning.  More  details 
will  be  discussed  in  subsection  (5). 

5)  Other  numerical  considerations 

The  forward  time  difference  and  center  space  difference  in  general 
are  used.  Except  for  , all  dependent  variables  are  evaluated  at  each  grid 
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point. Eq.  4 is  solved  l?y  the  block  relaxation  aethod.  Because  of  reflection 
of  spurious  energy  of  the  internal  gravity  wave  from  the  upjier  boundary 
originally  set  up  by  mountain  valley  terrain,  the  reflected  wave  energy 
leads  to  undesirable  computation  results,  ^he  spurious  wave  energy  can 
be  partly  eliminated  by  introducing  some  Rayleigh  friction  hear  the  upper 
boundary,  ^he  Rayleigh  friction  terms  may  be  written  as  -k(u-U),  and 
-kv  if  V = 0.  In  addition,  since  the  horizontal  eddy  exchange  terms  are 
not  included  in  the  model,  a simpler  smoothing  approach  is  used  instead. 

A three-point  smoothing  formula  for  each  variable  is  given  as  follows: 


where  s is  a constant  to  be  chosen. 

As  mentioned  in  the  preceeding  sub-section,  the  mountain  was  assumed 
to  be  flat  initially  and  will  be  gradually  built  up  over  26  time  steps 
until  the  final  desired  ridge  heights  are  reached.  The  undesired  initial 
disturbance  can  be  eliminated. 

III.  Computation  results 

Several  computer  runs  were  made  for  a combination  of  various  boundary 
conditions  for  both  wind  and  temperature,  different  geometric.-l  factors  of 
the  mountain-valley  terrain,  and  different  turbulence  paramatrization  in 
the  surface  layer.  Table  1 shows  the  various  basic  inputs  for  different 
runs.  In  each  computer  run,  the  isopleths  of  the  three  wind  components 
u,  V,  w,  and  T for  cross-sections  in  the  x-z  plane  at  J = 7 and/or  J = 4 
are  illustrated.  The  x-,  y-,  and  z-axes  directions  represent  the  cross- 
valley horizontal,  but  pointing  toward  upvalley,  and  vertically  upward 
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directions,  respectively.  To  be  explicit,  a positive  u is  the  horizontal 
wind  blowing  from  the  left  to  the  right  in  fig.  la;  a positive  v is  the 
horizontal  wind  blowing  upvalley  into  the  paper;  a positive  w is  the  vert- 
ical wind  blowing  bottom  upward.  The  cross-section  at  Jf  = 7 will  illustrate 
the  circulation  and  thermal  situation  for  the  middle  of  the  valley  and 
the  cross-section  at  J = 4 will  illustrate  the  information  at  the  lower 
part  of  the  valley.  At  J = 4,  the  valley  floor  and  side  ridges  start  to 
level  out.  Several  cases,  e.g.  Figs.  10,11,12,  show  that  the  height  of 
the  valley  side  ridges  above  the  valley  floor  start  to  drop  to  zero  at  J = 4, 
In  other  words,  there  is  a completely  flat  plain  between  J=1  and  4,  while 
a valley  slope  exists  from  J=4  to  10.  Between  J=10  and  11  it  is  a flat 
plateau.  The  horizontal  grid  size  and  height  of  ridge  above  the  valley 
floor  vary  fr9ra  case  to  case  although  the  same  diagram  is  used.  The  vertical 
grid  size  is  fixed.  Therefore  the  mountain  scale  is  somewhat  out  of  pro- 
portion with  the  mesh  scale  used  in  the  illustration.  The  isopleths  for 
u,v,w  are  isoteohs  of  the  corresponding  wind  components  in  m sec”^.  The 
isopleths  for  T are  temperature  deviations  from  horizontal  mean  in  °C. 

Since  the  gradients  of  these  quantities  can  be  very  large,  the  interval  of 
Isopleths  used  in  the  following  figure  are  not  always  uniform  and  all  iso- 
pleths will  be  labeled. 

Results  of  Run  # 1 are  shown  in  Fig.  3*  fhe  striking  feature  is 
the  tilt  of  positive  and  negative  axes  of  u,  w,  and  T toward  the  upstream 
direction  (negative  x-axis)  of  the  basic  flow.  The  maximum  +u  is  correlated 
very  well  with  the  -w  and  -T  6ver  the  lee  near  the  ridge.  In  3a  an  area 
of  negative  isotech  was  situated  near  the  center  of  the  valley  but  it  covers 
nearly  half  the  windward  slope.  This  signifies  the  development  of  the 
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return  flow  and  a separatedcell  near  the  bottom  of  the  valley,  when  the 
resultant  wind  In  the  x-z  plane  Is  computed  by  using  the  w-component. 

The  development  la  possible  because  a negative  and  a positive  w area  are 
situated  over  the  windward  and  small  part  of  lee  slope  In  the  bottom  center 
of  the  valley.  Owing  to  the  tilts  of  the  u and  w fields,  waves  in  the 
x-z  plane  at  J=4  and  ? with  tilts  axes  resulted,  similar  to  many  moxintain 
wave  patterns.  However  the  wavelength  is  larger  but  the  corresponding  wave 
amplitude  is  smaller  because  the  present  boundary  conditions  are  somewhat 
different  from  the  typical  investigation  of  mountain  wave  problems,  fer 
the  small  or  meso-scale  circulation,  the  periodic  condition  seems  to  be 
more  reasonable  than  the  boundary  conditions  used  in  the  mountain  irave  ’ 
studies  with  a single  mountain. 

Cross-soctions  of  the  flovre  in  the  valley  axis  direction  at  J=4  and 
J=7  are  shown  in  1 .b  and  1.f.  It  is  found  that  the  upvalley  flows  are 
maximum  near  the  center  of  the  valley  at  about  500m  above  the  valley  floor. 
In  comparing  the  magnitudes  of  v at  J=4  and  at  J=7,  the  magnitude  of  v 
component  is  greater  halfway  along  the  valley  (J=7)  than  near  the  valley 
mouth  (J=4),  Comparison  of  resultant  u and  w in  the  x-z  plane  show  there 
is  not  much  difference  in  strength. 

The  temperature  deviation  from  the  mean,  T,  for  cross-sections  at 
both  J=4  and  J=7  are  of  little  difference.  Again  the  inclination  is  to- 
ward the  upstream  direction.  Large  negative  T covers  most  of  the  \idndward 
slope,  whereas  a very  small  section  is  found  near  the  ridge  at  the  lee, 

Pbr  the  center  of  the  valley  from  the  ground  up  to  2km  are  positive  values. 
The  2-dimon3ional  results  obtained  here  are  reasonable  and  similar  to 
results  obtained  by  Tang  (1976). 
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Run  # 2 is  coaputed  for  a 24km  wide  valley  and  a valley  slope  of 
7x10“3  terrain  for  about  2 hours  and  is  shown  in  Pig.  4,  Except  for  the 
relative  narrowness  of  the  valley  and  the  steeper  slopes  both  in  the  valley 
sides  and  axis  directions,  the  assumptions  are  basically  the  same.  The 
patterns  of  u and  w are  similar.  The  wind  speeds  u at  lower  levels  are 
relatively  small.  The  v component,  however,  near  the  center  of  the  valley, 
is  about  0.30m  sec"^  higher.  The  level  of  maximum  v is  again  about  500m 
above  floor.  At  the  4 to  5 kw  level  there  is  a return  flow  in  the  down- 
valley  direction  and  it  appears  very  reasonable.  The  appearance  of  the 
T field  has  changed  somewhat  and  is  dominated  by  positive  deviations.  The 
negative  area  shrinks.  The  air  returns  warmer  in  the  valley  center  than 
over  valley  sides.  The  general  v/arming  is  believed  to  be  attributed  to  the 
stronger  downward  motion  above  the  steeper  slope  of  the  ground . 

An  even  narrower  and  steeper  valley  is  used  as  a terrain  model  for 
the  study.  The  width  is  4 km  and  the  valley  slope  is  33x10"^.  dines  the 
time  interval  used  for  the  integration  cannot  be  larger  than  5 sec,  only 
about  29  minutes  (350  steps)  of  computation  have  been  made  for  this  run. 

The  appearance  of  the  patterns  have  changed  somewhat.  The  axes  become  loss 
Inclined  than  the  first  two  cases.  The  u component  near  the  center  of  the 
valley  becomes  even  smaller  than  that  near  the  ridge  level.  The  return 
flow  region  remains  in  the  center  of  the  valley  and  over  the  windward  slope. 
The  positive  and  negative  w fields  become  distributed  more  or  less  verti- 
cally over  the  windward  and  leeward  slopes  respectively.  The  center  of 
maximum  magnitude  is  about  2km  above  both  slopes  between  the  ridge  and  the 
center  of  the  slopes.  The  maximum  magnitude  is  about  3*5  >*>  aec"^.  There- 
fore the  troughs  or  the  wave  in  the  x-a  plane  will  bo  more  or  less  vertical 
The  maximum  amplitude  of  the  wave  is  about  2km  above  valley  side  slopes. 
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The  T field  has  large  negative  values  near  the  center  but  is  warmer  than 
over  the  slopes.  The  coldest  spot  in  the  air  is  near  the  top  of  the  wind- 
ward slope.  The  flow  in  the  horizontal  direction  of  the  valley  axis 
direction  is  largest  for  the  last  eases.  The  maximum  valley  wind  is  about 


1,5  n sec”^  between  20Cm  and  500m  above  the  ground,  the  maximum  down-valley 
wind  is  about  1.4  m sec~^  at  about  5 lo  ^ loa  above  ground. 

The  first  2 cases  are  computed  for  a sloping  valley.  The  top  and  bottom 
of  the  valley  are  assumed  to  bo  horizontal  by  assuming  that  the  variations 
of  all  quantities  in  the  y-axis  direction  are  zero,  fbr  the  following  cases 
the  boundary  conditions  in  the  y-axis  is  again  assumed  to  bo  cyclic.  Fbr 
Runs  4,  5.  6 and  7 the  height  of  the  ridges  above  the  valley  floor  is  cons- 
tant, but  it  may  be  a different  constant  for  different  cases.  Other  runs 
are  for  cases  where  the  height  of  the  valley  side  ridges  start  to  drop  off 
from  position  at  J=7  to  zero  at  J=4,  a typical  valley- plain  topography. 

Runs  #4  and  #5  are  for  valley  with  a 64  km  wide  valley.  Run  |4  is  made 
with  a warmer  valley  slope  which  cools  with  an  assumed  rate  as  time  goes 
on.  Results  are  illustrated  for  cross-section  at  the  middle  of  the  valley 
slope  where  J=7  at  time  20100  sec  ( 5 hrs).  Run  #5  is  made  with  a heating 

case.  The  axes  of  these  patterns  are  slightly  inclined  downstream  of  the 
basic  flow.  However,  since  the  vertical  velocity  field  w is  very  weak 
and  is  essentially  vertical,  the  wave  trough  will  remain  vertical, and 
energy  propagating  downward  is  negligible.  The  characteristics  of  the 
return  flow  for  those  cases  are  again  in  agreement  with  earlier  findings. 

The  intensity  under  the  cooling  condition  is  stronger  than  that  under  the 
heating  condition.  The  flow  along  the  valley  is  toward  the  down  and  up 


valley  directions  for  Runs  #4  and  #5  respectively.  This  is  reasonable. 
The  down-valley  component  (-3*5  seo“^)  in  the  cooling  case  is  stronger 
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than  the  up-valley  flow  (+2.4  m sec”^).  This  is  in  general  agreement  with 
observational  results  obtained  in  Vermont  valleys* 

Similar  to  the  assumption  used  in  Runs  4 and  5,  Runs  6 and  7 are  com- 
puted with  a valley  only  half  as  wide  as  the  valleys  in  Runs  4 and  5*  In 
comparing  these  two  sets,  again  patters  are  very  similar  to  corresponding 
eases  and  a similar  conclusion  can  also  be  drawn.  Between  Runs  4 and  5 the 
cooling  surface  condition  leads  again  to  stronger  intensity  than  with  the 
heating  condition. 

Runs  #8  and  #9  are  computed  for  u,  v,  w,  and  T at  position  J=7  under 
the  same  assumption  except  for  the  corresponding  surface  cooling  and  heat- 
ing conditions  on  a realistic  valley-plain  terrain.  The  variation  with 
time  in  each  case  are  illustrated  by  the  results  obtained  at  t=10100  sec 
and  at  t=20100  sec.  The  one  with  the  cooling  surface  condition  produces 
stronger  intensity  than  the  one  with  the  heating  surface  condition  in 
general.  As  time  increrses  to  about  5»5  hrs,  nothings  changes  noticeably. 
The  only  quantity  that  is  intensified  is  the  upvalley  or  downvalley  flow. 

Run  #10  is  the  experiment  with  surface  layer  formulation  based  on  the 
equations  used  in  Section  II. 3.  Since  the  time  step  has  to  be  very  small 
when  the  surface  layer  formulation  is  included,  computation  has  been  made 
for  750  sac.  Results  are  shown  in  i'ig.  12.  The  patterns  of  the  isopleths 
are  very  similar  to  those  shown  in  Fig.  4 except  that  the  basic  wind 
profile  in  the  vertical  direction  decreases  with  height  in  the  upper  layer, 
which  can  be  attributed  to  the  input  profile  and  small  mountain  slope  used 
in  Fig.  4.  It  appears  that  using  the  surface  layer  formulation  has  made 
little  difference  upon  the  gross  features  of  the  circulation.  However, 
further  experiments  including  the  surface  layer  should  be  pursued. 
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The  present  numerical  model  has  produced  several  Interesting  and  en- 
couraging results  of  a three-dimensional  model  for  mountain-valley  vd.nd 
systems.  E(y  no  means  is  tnis  an  ultimate  model  for  prediction.  While 
pursxiing  the  Improvement,  some  results  obtained  here  may  be  used  for  com- 
parison with  observations  available  at  hand.  Based  on  studies  by  Davidson 
(i960),  for  a valley  with  ridge  height  of  800m  the  maximum  winds  are  found 
at  about  400m  level  which  is  halfway  to  the  ridge  level.  The  comp'rison 
of  tne  level  of  maximum  down-valley  wind  with  the  results  obtained  here  is 
shown  in  fig.  13«  lAth  the  assumed  topography  used  here  we  have  obtained 
the  maximum  down-valley  wind,  generally  occuring  at  the  5^0®  level  which 
is  in  close  agreement  with  observation.  Even  for  many  computations  the 
down-valley  wind  speed  is  small,  . the  level  of  maximum  remains  in  the 
neighborhood  of  the  middle  level  of  the  valley.  In  the  slope  wind  theory 
for  an  infinitely  long  slope  the  level  of  maximum  wind  is  only  about  50m 
above  the  ground.  It  cannot-  be  used  to  explain  the  valley  wind  phenomenon. 

A simple  three-dimensional  linear  theory  cannot  explain  this  phenomenon  either 
It  is  believed  that  the  non-linear  theory  has  a contribution  to  this  con- 
clusion. 

IV.  Discussion 

The  numerical  model  developed  in  this  study  has  included  the  non- 
nydrostatic  assumption  and  turbulence  parameterizations  that  are  believed 
to  be  important  in  the  study  of  mountain-valley  wind  systems  of  steeper 
terrain.  The  computations  yield  reasonable  results  in  the  three-dimen- 
sional flow  over  a typical  mountain- valley  terrain.  The  mountain-valley 
wind  system  resulted  from  the  non-linear  interaction  between  the  thermally 
Induced  flow  in  the  valley  and  the  upper  air  flow  of  the  larger  scale  motion. 


This  model  produce^  the  following  conclusions: 

1)  The  positions  of  the  return  flow  or  separated  cell  in  the  valley  obtained 
in  this  study  is  in  agreement  with  the  earlier  study. 

2)  The  imposed  nighttime  thermal  surface  condition  produces  down-valley 
wind,  the  imposed  daytime  thermal  surface  condition  produces  up-valley  wind. 

3)  The  maximum  dovm-valley  wind  level  is  about  500m,  for  a ridge  1km  above 
the  valley  floor. 

4)  The  up-valley  vdnd  is  relatively  weak  as  compared  to  the  down-valley 
wind  speed  while  the  other  conditions  remain  the  same. 

5)  The  component  of  cross-valley  wind  is  maximum  at  the  lee  near  the  ridge, 
or  above  the  ridge. 

6)  The  wave  flow  pattern  in  the  upper  part  of  the  troposphere  is  similar 
to  some  earlier  findings . 

In  the  future  studies,  continuous  effort  will  be  placed  on  the  improve- 
ment of  the  accuracy  of  the  model.  This  m?y  be  corrected  by  modiiying 
the  present  computation  scheme  and  increasing  the  numbers  of  the  horizontal 
grids.  In  order  to  be  more  realistic,  cyclic  boundary  conditions  will  be 
modified  so  that  they  can  be  applied  to  actual  conditions,  first  we  will 
conduct  a study  based  on  an  analytical  model  of  two  mountain  ridges  and 
sloping  valleys,  then  later  compare  the  results  from  the  numerical  and 
theoretical  models.  Actual  observational  data  from  the  Tularosa  Basin  will  be  used  to 
validate  the  results.  Finally  an  attempt  will  be  made  to  build  a prediction 
model  for  the  1\hite  Sands  Missile  Range. 
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Table  1 


Basic  inputs  for  different  computer  runs 


The  number  Maximum  ridge  height 
of  run  and  above  valley  floor 
time  (sec)  and  floor  height 
above  plain  (m) 


Geometrical  factors  of  terrain 


Surface  thermal  condition 
Cooling  I Neutral!  Hea  ting 


o 

o 

1000+300 

2 (7550) 

1000-300 

3 (1750) 

1000+6:0 

4 (20100) 

ICOO+6OO 

5 (20100) 

1000+600 

6 (22500) 

1000+300 

7 (22500) 

1000+300 

8 (10100- 
20100) 

500+150 

9 (10100- 
20100) 

500+150 

10  (750) 

1000+600 

w 


rldg« 


A*or  ^ 

■ 1 


T 


1 1 


3 5 7 S ? 


ng.  2.  The  height  of  ridge  and  height  of  the  valley  floor  in  the 
y-z  plane  at  the  middle  of  the  valley. 
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1 . Introduction 


The  diurnal  surfaco  thermal  condition  has  been  shovn  to 
bo  an  important  controlling  factor  on  daily  circulation  in  the  pla- 
netary boundary  layer. 

Recent  investigations  of  mountain-valley  circulation  by 
Tang(1976),  and  by  Tang  and  Peng  (197^)  with  a numerical  model,  show 
that  thermal  boundary  conditions  have  an  essential  influence  on  cir- 
culation in  the  boundary  layer.  The  assumption  on  the  thermal  bound- 
ary condition  for  their  investigation  is  a prescribed  surfaco  temp- 
erature or  heat  fl\ax  at  the  ground  surfaco.  Since  the  moiintain-valley 
terrain  consists  of  \mevon  surfaces  such  as  slopes,  ridges,  and  valleys, 
the  amount  of  solar  and  long-wave  radiation  over  those  surfaces  is 
different  from  that  on  a horizontal  plane.  Incorporation  of  a variable 
surfaco  heat  fliox  in  the  study  would  make  it  a very  complex  problem.  A 
most  realistic  mountain-valley  circulation  model  considering  this  temp- 
orally and  spatially  variable  surface  heat  flux  is  yet  to  bo  developed. 
The  purpose  of  this  study  is  to  derive  an  analytical  expression  of  the 
ground  surfaco  temperature  and  heat  flux,  at  each  small  time  stop,  in 
terms  of  solar  and  long-vrave  radiation,  geometrical  factors  of  the  ter- 
rain, and  initial  surfaco  temperature.  The  result  can  be  utilized  for 

4 

a numerical  experiment  of  atmospheric  motion  and  temperature  over  a 
sinpl.e  mountain-valley  terrain. 

In  two  previous  studies  Rah  (1959,  19^^)  invesitgated  the 
influence  of  a sloping  surface  on  radiation,  the  present  study  develops 
a method  for  obtaining  surface  temperature  by  solving  the  equation  of 


1.  Introduction 

The  diurnal  surfaco  thermal  condition  has  been  shovjn  to 
be  an  important  controlling  factor  on  daily  circulation  in  the  pla/- 
netary  boundary  layer. 

Recent  investigations  of  mountain-valley  circulation  by 
Tang(1976),  and  by  Tang  and  Peng  (197^)  with  a numerical  model,  show 
that  thermal  boundary  conditions  have  an  essential  influence  on  cir- 
culation in  the  boundary  layer.  The  assumption  on  the  thermal  bound- 
ary condition  for  their  investigation  is  a prescribed  surfaco  temp- 
erature or  heat  flux  at  the  ground  surface.  Since  the  movuitain- valley 
terrain  consists  of  rxnoven  surfaces  such  as  slopes,  ridges,  and 'valleys, 
the  amount  of  solar  and  long-wave  radiation  over  these  surfaces  is 
different  from  that  on  a horizontal  plane.  Incorporation  of  a variable 
svurfaco  heat  flux  in  the  study  would  make  it  a very  complex  problem.  A 
most  realistic  mountain-valley  circulation  model  considering  this  temp- 
orally and  spatially  variable  surface  heat  flux  is  yet  to  bo  developed. 

The  purpose  of  this  study  is  to  derive  an  analytical  expression  of  the 
ground  surface  temperature  and  heat  flux,  at  each  small  time  step,  in 
terms  of  solar  and  long-vrave  radiation,  geometrical  factors  of  the  ter- 
rain, and  initial  surface  temperature.  The  result  can  be  utilized  for 
a numerical  experiment  of  atmospheric  motion  and  temperature  over  a 
simple  mountain-valley  terrain. 

In  two  previous  studies  FUh  (1959.  19^^)  invesitgated  the 
influence  of  a sloping  surface  on  radiation,  the  present  study  develops 
a method  for  obtaining  surface  temperature  by  solving  the  equation  of 
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heat  conduction  of  the  soil  and  the  energy  balance  equation  at  the  ground, 
•by  including  the  terrain  influence  on  radiation  . The  heat  balance 
condition  at  the  ground  surface  over  a mountain-valley  terrain  at  dif- 
ferent locations  across  the  valley  is  conceivably  not  the  same.  As  a 
consequence  the  ground  surface  temperature  will  vary  from  place  to  place, 
at  least  in  the  daytime.  The  heat  flux  from  the  ground  surface  id.ll 
Interact  vdth  atmospheric  motion  and  temperature  to  develop  a very 
complicated  dynamic  system.  Our  ultimate  goal  is  to  solve  this  com- 
plete system  numerically  for  mountain-valley  circulation,  Hoi<raver,  our 
immediate  interest  is  the  formulation  of  the  ground  surface  temperature 
of  a typical  valley  terrain, 

2.  Solar  radiation  on  the  sloping  surface 

Pbr  an  undulating  mountain- valley  terrain  as  shoim  in  Rg.  1, 
the  mountain  ridge  on  the  east  side  of  a north-south  valley  idll  cast 
its  shadow  on  the  west  side  of  the  valley  (Fig.  1).  The  solar  radia- 
tion reached  on  the  slope  may  be  written  as 

S = Sq  cos  i , (1) 

where  Sq  = solar  constant 

i = the  acute  angle  between  the  sun's  rays  and  the  normal  to 
the  slope  n , 

Then  cos  i can  be  written  as 

COS  I = cos  Of  sin  h + slnacos  hcos  ^ ^ (2) 


where  a = the  angle  of  inclination  of  the  slope 
h = the  elevation  angle 


♦ = Vo“''n 

4>o, 


if  measured  counter-clockvdsc  from  the  south; 
if  measured  clockv.'ise  from  the  south) 


2 


WEST  < > east 


Fig  . t , A projective  geometrical  relationship  between  the' sun’s 
ray  and  the  terrain.  • ■ - 
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A s azimuth  of  the  sun 

♦jj  = azimuth  of  the  normal  to  the  slope. 

Using  the  laws  of  cosine  and  sine  in  spherical  trignometiy  and  the 
geometrical  relationships  of  various  angles  shown  in  ELg.  2,  one  ob- 
tains the  follovdng  set  of  equations  (e.g.,  see  Robinson,  I966) 


felnh  = sln9sln6  + cos  7 cos  6 cos  cu 

coshsini^Q  = cos  6 sin u) 

cosh  cos  9 cos  = sin  (p  sin  h - sin  6 


(3) 

W 

(5) 


where  6 = latitude 

9 • s:  declination  of  the  sun 
m = hour  angle. 

Substituting  (3),(^),  and  (5)  into  (2),  eliminating  and  rear- 
ranging yields  the  radiation  flux  on  the  slope. 

S “ Jsinfi  + K cos  6 cosu)  + sini^jjSino  cos  6 sinu))  (6) 


where  *j  = sin9COsa.  - cos  9 s in  a cos 
IC  cos  9 cos  or  + sin9sinorcosijrjj 


(7) 


fbr  an  east  (or  vrest)  -facing  slope,,  ~ 

Sg  ^ (sin9cos«sin6  + cos  9 cos  cr  cos  i cos  w 


IT 


4-  slnacbs  6 sintu) 


fbr  a horizontal  surface,  or  = 0 


*=  S = S. (sin 9 sin  & + cos  9 cos  6 cos  uu) 


S — W w* 

e,w  ® 


(8) 


(9) 


This  can  be  applied  to  the  level  surfaces  on  the  ridge  and  valley  floor, 
fbr  a case  on  the  valley  floor,  at  point  A in  H.g.  1 for  example,  a 
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critical  hour-angle  a>  corresponding  to  sunrise  or  sunset  has  to  be 

c 

determined. 

Since 


cos  t 


V 


X * + V * 

yg 


Xg*  + Zg* 


then  sin  = 


+ yg* 


cos  ^ 

■ ■ 

cot*hg^  + cos^ij 


Ibr  the  present  case,  •where  the  ■valley  axis  is  in  the  north-south  di- 
rection, f n *s, 

Z * 

Equation  (12)  then  becomes 


Sin  hg  = 


^cot*ti_.  + sln*i 


Since  ♦0''  cannot  be  determined  tdien  the  corresponding  o)  is  not  known, 
one  cannot  solve  (13) • Therefore  we  need  another  equation  to  determine 
the  relationship  between  to  and  • This  can  be  achieved  by  using  the 
relationships  (3),  (4),  and  (5)  to  obtain 


8ln*t 


sin^o) 


sin^io  + (cos  tp  tan  8 + sin  9 cos  uj)* 


The  critical  u>^  can  bo  obtained  by  the  ’cut  and  try’  method  vdth  given 
9 and  t from 
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sln*a) 


sin^u)  + (cos  cp  tan  6 + sin9cosu))’ 


cot*h + 


sln*<u 


SPC  sln^tu  + (cos  9 tan  6 + sin  9 cos  uj)^ 


|sln9 


sln&  + cos  9 cds  6 cos  ( 


Another  method  to  determine  (u^  is  to  simply  set  Eq.  (8)  to  zero.  How- 
ever, the  angle  or  has  to  be  changed  into  the  appropriate  critical  so- 
lar elevation  angle  for  the  point  in  question. 

^y  this  method  the  solar  flux  of  all  positions  on  the  mountain- 
vall^  surface  can  be  computed.  If  tho  physical  surface  of  the  motin- 
tain-valley  terrain  is  symmetrical  and  periodic  in  space,  the  total 
energy  received  in  24  hours  will  be  the  same  for  both  east  and  west 
slopes.  However,  the  accumulated  energy  received  at  symmetrical  positions 
on  opposite  slopes  will  be  different  at  any  given  instant  during  the  day. 
Jbr  practical  purposes  one  may  just  compute  the  flux  every  half  hovu* 
for  some  typical  locations,  fbr  further  computations,  the  data  can  be 
Interpolated  for  any  time  interval  and  position  needed.  This  flux  id.!! 
be  used  to  determine  the  surface  temperature  and  heat  flux. 

3.  Long- wave  radiation 

The  next  most  important  item  in  the  energy  balance  equation 
after  solar  radiation  is  the  effective  long-wave  radiation  flux.  Because 
the  mountain  surface  itself  can  screen  away  the  IR  radiation  for  a 
given  azimuth  ^ , and  a given  elevation  angle  h subtended  by  the  mountain, 
tho  effective  IR  flux  can  be  written  as  (Rih,  1964): 


I 


Jt 


w 

rTT 


T • ^ ^ - e ) cos  I cos  h 

* “ s g 


hdh 


W) 


.m 


A--VT  # 

+ \ d*  \ (fis  - e_)  cos  l.cos 

\c») 


h dh 


(16) 


vhere  eg  = IR  radiation  intensity  on  the  slope  at  the  point  in  question 
eg  = IR  radiation  from  the  atmosphere 
Og  = IR  radiation  from  the  opposite  slope 
V*)  = elevation  angle  produced  by  the  slope  surface 
(inclination  ® ) itself  in  the  direction  . 


In  this  study,  ridge  lines  are  assumed  to  be  parallel  as 
shovn  schematically  in  figure  1 , and  h(  ^ ) for  the  elevation  angle 
of  east  and  vest  slopes  observed  at  Point  A,  hQ  and  hg  respectively, 
can  be  vritten  as  . 


arc sin 


cos 


2 4, 


cot2(hg^)  + cos2  ^ 


and 


hg  s arcsin  ^ 


cos2  t|» 


A 

2 


cot2(hgjj)  + cos2  »[» 
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The  anglo  ( <•  ) is  tha  elevation  angle  due  to  the  slopo  itself  where 
the  observation  site  stands. 

If  the  IR  radiation  is  isotropic  xdthin  the  integration  Units 
and  the  valley  is  infinitely  long,  the  integral  in  Eiiuation  (16)  can  be 
evaluated  and  becomes 


vdiere 


*=  .1  ( Co  - e ) S + * ( e„  - ) ( 1 - ^ S ) 


( (cos(hgx)+cos(h^))*cosa  - (5ln(hgx)-5in(h-J.sin  a } 


bben  the  black  body  tosiperaturos  are  used,  I at  A is 


• cos  a - (sin  h - sin  h^v)  sin  a + ct(0  X*  - 3 ) (21) 

I s s e e ' • 

where  Ijj  s effective  IR  radiation  flux  (i.e.  difference  of  IR  radiation 
between  mountain  surfaces  and  atmosphere)  on  a horizontal 
surface 

T5  = the  absolute  temperature  of  the  surface  at  the  point  in 
question 

Tjj  » the  absolute  temperature  of  the  opposite  slope 
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^ *=  the  enisslvlty  of  the  earth  surface  at  the  point  in 
question,  and  of  the  opposite  slope,  respectiveJy 
hgx  . l*cx  ® ®lovation  angle  due  to  mountain  ridges  GH  and  CD  respectively 
C s Stefan-Boltzmann  constant. 

Since  or  s 25  as  assumed  in  the  present  model  and  hg^  £ a » and  , 

it  is  evident  that 

sin  eif(sln  - sin  « (cos  hgj^  + cos  cos  of ; 

furthermore  T^-T^  is  gensrally  small.  The  overall  affective  IR  radia- 
tion may  bo  apprbximated  as 

Si  V ^'cx^  • (22)  . 

However  the  quantity  , defined  as  the  effective  IR  radiation  on  a 
horizontal  surface,  must  be  calculated  iteratively.  We  nay  use  an 
euqjirical  formula  for  Ijj  developed  by  Brunt  (1939),  i.o., 

Ijj  = cTg^  (1-a-b^)^  (23) 

idiere  a = 0.^ 

b a 0.08 

o = surface  vapor  pressure  in  millibars 
^SA  s=  the  surface  temperature  at  the  point  in  question, 

(In  the  present  problem  Tg  is  the  temperature  at  the  point  A,  which 
is  used  to  calculate  Ij^  at  A),  finally  we  may  write  the  effective  IR 
radiative  flux  as 

I « -^oTg  (1  - a - b^}(cos  hg^  + cos  cos  or.  (24) 

The  value  of  I depends  on  the  position  on  the  mountain-valley  surface 
and  can  be  calculated  for  each  point  on  the  surface.  The  surface  temp- 


10 


I 


erature  in  the  expression  for  I remains  to  be  determined  later  by 
Bolvlng  the  heat  conduction  equation  of  the  soil  and  a system  of  thezmo- 
lydrodynamic  equations  of  the  atmosphere  simultaneously. 

Heat  conduction  equation  of  soil  and  the  boundary  conditions 

• The  general  equation  of  the  mountain- valley  surface  may  be  writ- 
ten as 

D D(x) 

Using  z to  indicate  the  height  above  the  fixed  level  surface,  and 
the  height  above  the  ground  surface,  we  have 
• z = D(x)  + C 

The  heat  conduction  equation  which  governs  the  variation  of  the  soil 
temperature  may  be  written  approximately  as 

’ET  “ 

where 

= thermal  diffusivity  of  the  soil 
Pg  = density  of  the  soil 
Og  s specific  heat  of  the  soil 

The  boundary  condition  at  the  ground  surface,  after  neglecting 
seatteidng,  precipitation  and  evaporation,  and  assuming  the  normal  temp- 
erature gradient  to  be  about  the  same  as  the  vertical  temperature  gra- 
dient, will  be  approximately  - - 

S - I + X -g-  - X,  IP  = 0 ^ (25) 

where  S = solar  radiation  reached  on  the  slope 

I - long-v;ave  exchange  between  tho  slope  and  the  environment 
* Xg=  thermal  diffusivity  of  tho  soil 
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X 


« eddy  conductivity  of  the  air  (assumed  constant  in  the 
layer  T.4m  above  ground,  and  may  vary  with  time) 

I , Tj  « temperatures  of  the  air  and  soil  respectively. 

Jbr  the  boundary  condition  at  Q » -•  , Tg  is  finite , 

5*  Solution  of  surface  temperature 

The  differential  equation  of  the  heat  conduction  of  the  soil 
■ay  be  written  in  finite  difference  form  for  time  differential.  Assume 


T,  - % + Bz  + T^.  (c,x,t)  , 


(26) 


where  B = prescribed  for  2>0, 
and  ' B = 0 for  z<0. 

Vb  designate  the  elevation  of  the  valley  floor  as  2=0,  Substituting 
Tg  into  the  finite  difference  equations  of  heat  conduction  of  soil  leads  to 


m + At  _ ,t 
At  • Ps 


(27) 


in  which  z is  replaced  by  C near  the  surface.  The  solution  of  this 
equation,  obtained  ly  varying  the  parameters  and  assuming  that  Tg*^  is 
a known  function,  and  can  be  written  as 


T/'*‘^^(x.C,t)  = A(x)e^^  + Be“^^ 


+ ■ v(C  - C’)  (-Y*T^*)dC*  (28) 


where  ^ 


'-^At  . X /c  p 5 • Applying  the  boundary  conditions  at  z - 

leads  to  B'  = 0.  Next  the  coefficient  A(x)  should  be  determined  from 


the  energy  balance  equation  at  C = ® with  analytical  solutions  of  S 
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*nd  I at  each  grid  point. 

Assuming  tho  condition  of  tha  soil  tomporatur©  at  the  air-soil 
interface  to  be  continuous,  we  have  the  detailed  expression  of  the 
energy  balance  equatioh  as  ' . 


t+At  „t+At 
T - Tg 


^ (sin  i?  cos  a - 


cos  tp  sin  or  cos  4 ) 


• sin  6 + (cos  (p  cos  or  + sin  «p  sin  a cos  • 


,t+At 


• cos  6 cos  0)  + sin  ^j^slnacosfi  slna> 


(cos  + cosh  ) cos  or 
' gx  cx 


Substituting  % - T„  + Bs  + n(x,s,t),  (25)  ^ rearranging 


leads  to 


+1^  + 2 aTo^  (1  - a-  b^e)  (cos  + cos  h^)cos 


idiere 


F(a.»„,a.)'  = _ x^B+  ^ (T,  + Br(x)) 

+ So  ((sin  9COSer  - cos  <p  slnof  cos  sin  6 + 

****** 

+ (cos  9 cos  or  + sln9  ^In  a cos^r^^)  cos  a cos  a> 
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The  soil  heat  JClux  at  the  ground  surface 


After  substituting 


0 and  solving  for 


2 or  Tq’  (1  -a  - b->/e)  (cos  h + cos  h ) cos 

LiX 


Ibr  a given  location  x and  time  t,  A is  constant  and  + At  ^ 

is  determined.  V7ith  known  T T can  be  calculated  if  the  t) 


erature  of  the  air  above  the  surface  is  kno^-m.  Set 


fbr  integration,  the  time  interval  t=  30  minutes.  Kg  = .75x  10"^  cra^sec'^ 
and  £ 0.86  cm“\  Diurnal  variation  in  soil  temperature  at  a depth 
of  50  cm  is  nearly  zero  at  O’Neill,  Nebraska,  Aug,  31,  1953,  based  on 
data  obtained  by  Lettau  and  Davidson  (1957).  It  follows  that  1 ^ 

50 

This  means  that  1 is  negligible  in  comparison  to  V , Neglecting  1 
all  together  from  (3^)  yields  an  error  loss  than  3/^.  meso-scalo 

calculations  this  error  is  acceptable.  Then  Tg’  may  be  simply 

written  as 


.t+6t 


‘s 


,t+At 


C=0  . XgV  + L 


+ 


- • (T  ' ) 

XsY  + L ^ s (36) 


KLnally  ve  have 

.t+  ht 

‘C=0 


= To  + Bz  + 


• pt+&t 


4.  - •(■p 


XsY+  L XgV  + L s 


Eq,  (36)  should  be  solved  simultaneously  vjith  the  set  of  finite 
difference  equations  governing  the  air  temperature  and  motion  in  the 
boundary  layer.  The  solution  can  be  obtained  if  the  air  temperature 
at  a height  az  above  the  surface  at  time  t +At,  and  the  ground  sur- 
face temperature  at  time  t,  are  knovm.  An  initial  rough  but  realis- 
tic estimate  of  air  temperature  + ^t  made  for  the  first 

computation  of  Tg^  +^t^  Vfi.th  the  flux  at  time  t +£*t  one  can  com- 
pute the  atmospheric  temperature  +2/^t  from  ^j^e  set  of  governing 
differential  equations  in  the  atmosphere.  Continue  this  computation 
iteratively  to  obtain  where  n is  an  integer  one  wishes  to 

assiune.  This  approximated  method  was  developed  mainly  for  the  purpose 
of  mountain- valley  circulation  studies.  It  can  also  be  applied  to 
studies  of  micrometeorological  problems  as  long  as  the  geometry  of  the 
terrain  is  simple.  In  the  present  model  the  point  of  interest,  the 
effective  IR  radiation  on  the  surface  Ig,  is  not  directly  calculated, 
but  rather  is  replaced  by  an  empirical  formula.  So  long  as  the  sur- 
face geometry  is  simple,  a simple  anlytical  solution  can  be  derived. 


' 5>  Results  and  Discussion 

Equation  (36)  Is  U3<jd  to  compute  ground  tomporaturo  for  three 
different  positions  In  a valley  located  5 mllos  north  of  Manchester  Cen- 
ter, Vermont,  They  are  Point  A on  tho  valley  floor.  Point  B on  the  west 
slope,  and  Point  F on  tho  east  slope.  The  heights  of  Points  B and  F az>e 
about  440m  arjd  420m  above  Point  A respectlvoly.  Tho  ridge  linos  are  548m 
above  the  valloy  floor  and  are  7»8kra  apart.  Point  A Is  approxlmatoly  2kin 
horizontal  distance  from  tho  west  ridge.  Tho  average  slope  at  Point  F 
is  rather  gontlo,  bolng  6.3°.  The  slope  on  tho  west  side  varies  from  12,5® 
on  tho  lou-or  half  to  18.8®  on  tho  upper  half  whore  Point  B Is, 

Other  physical  parameters  used  in  tho  computation  are  as  follows; 

X c 0.239  cal  cnr^soc“^dog“^  a =«  0.44 

X = also  varies  as  a function  of  b = O.O8 

time(soo  Table  1) 

XfiSt  0,001  e B to 

t - 1800  sec  z B 140  cm 

Cg  Pg  =0,4  cal  cm"^dog"^  6 = 20® 

albedo  = 0.15  ^ 9 = 43,16® 

Table  1 . Assumed  variation  of  \ with  time  over  tho  course  of  a day. 

hour  0 2 4 6 8 10  12  14  16  18  20  22  24 

(EST) 

X (cal  cm" ^ .05  .05  .06  .10  .24  .28  .30  .30  .20  .12  .08  .05  .05 

deg-^soo"^) 

The  sunrise  and  sunset  hours  for  Points  A,B,  and  F computed 
from  this  model  are  (0515.  1752  E.S.T.),  (0445.  I735  E.S.T.),  and  (05I5, 
1911  E.S.T.)  rospoctivoly. 
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The  air  temperatures  at  the  three  points  are  recorded  by  theino- 
graphs  located  1 «4ffl  above  the  gound.  The  data  of  August  7,  1957  taken 
f^a  Davidson  and  Rao  (1953)  used  in  the  computation  and  are  shown 
In  figure  3 dashed  lines. 

Since  there  is  no  observed  groxmd  temperature  Tg  at  each  site, 
a reasonable  initial  Tg  or  T^'  value  has  to  be  given  for  input.  This 
value  is  determined  with  a trial  run  using  an  appropriate  Tg  or  Tg'  un- 
til the  result  converges  to  a constant  and  its  rate  of  change  is  about 
the  same  or  slightly  larger  than  that  of  air  temperature  at  midnight  for 
a given  set  of  assumed  values  for  parameters,  fbr  daily  use,  Tg  can  be 
observed  and  there  is  no  need  for  such  trial  runs. 

Based  on  the  physical  configuration  of  the  terrain  and  input 
air  temperatures  for  each  point,  the  corresponding  computed  ground  temp- 
eratures and  upward  heat  fluxes  under  three  different  assumptions  of  A 
are  shown  in  Figures  3 and  4 respectively.  In  Figure  3 the  computed  ground 
temperatures  for  cases  of  constant  A , variable  X , and  doubled  value  of 
the  va Triable  A are  denoted  by  light,  heavy,  and  dash-dotted  curves.  The 
observed  air  temperatures  are  denoted  by  dashed  curves.  In  Figure  4 the 
computed  upward  heat  fluxes  for  constant  A , smaller  variable!  , and 
larger  variable  A cases  are  denoted  by  dashed,  solid,  and  dash-dotted 
curves  respectively. 

The  computed  sunrise  and  sunset  hours  shown  earlier  are  correct 
because  two  different  formulas  mentioned  earlier  give  the  same  computed 
results.  Besides  the  computed  sunrise  hour  at  Point  B is  the  earliest, 
and  the  computed  sunset  hour  at  Point  F is  the  latest,  among  the  three 
sites.  All  of  those  results  are  in  close  agreement  with  values  from  other 
sources. 
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figure  3.  Obser\red  diurnal  air  teaporature  (dashed  curve), end  coaputod  diurnal 
ground  toaperature  for  variable  \ (heavy  solid  curve),  for  constant  X (light  solid 
curve),  and  for  doubled  variable  X (dash-dotted  curve),  for  Point  A on  the  valley 
floor.  Point  B on  the  vast  slope,  and  and  Point  F on  the  east  slope  of  the  valley. 


Examination  shows  that  the  computed  groxmd  temperature  reaches 
the  level  of  air  temperature  at  each  site  about  half  an  hour  after  sun- 
rise. Then  both  air  and  ground  temperatures  start  to  rise.  The  times 
before  attaining  maximum  T^  of  the  Individual  site  occur  first  on  the 
west  slope,  then  on  the  valley  floor,  and  finally  on  the  east  slope  in 
figure  3»  This  is  as  expected  because  of  the  present  topographic  config- 
uration. Tg  starts  to  drop  generally  in  the  afternoon  but  drops  preci- 
pitously after  sunset  for  Point  A and  just  before  sunset  for  Points  B 
and  F.  The  phenomenon  at  Points  B and  F can  probably  be  attributed  to 
the  large  drop  in  air  temperature,  and  to  turbulence  taking  aaray  heat 
from  the  ground  and  reducing  ground  temperature  much  faster.  Another  in- 
teresting phenomenon  apparent  from  figure  3 is  that  the  daily  range  of 
ground  temperature  is  largest  at  the  site  in  the  valley,  being  a few  de- 
grees greater  than  at  the  sites  on  the  slopes.  The  differences  la  the 
ranges  of  ground  temperature  are  in  agreement  with  those  of  air  tempera- 
ture. This  phenomenon  appears  in  many  valleys. 

The  results  of  the  heat  flux  computations  in  Figure  k have  al- 
most the  same  shape  and  magnitude  for  the  three  different  cases  of  X . 
Although  the  ground  temperatures  Tg  for  different  X 's  are  different  at 
a particular  site,  the  changes  in  heat  flux  are  slight.  The  reason  is 
that  the  product  of  the  smaller  x with  the  large  vertical  temperature 
gradient  (between  the  air  temperature  at  1 ,4m  level  and  ground  temperature) 
Is  nearly  the  same  as  the  product  of  the  larger  x with  the  smaller  vert- 
ical gradient.  The  system  adjusts  itself  so  that  the  change  in  flux  is  gra- 
dual. The  trend  in  the  order  of  occurence  of  maximum  heat  flux  follows  that 
of  the  order  of  occurence  of  maximvim  temperature  at  tho  three  sites.  The 
total  heat  flux  wiiich  diffused  due  to  edd^'  processes  was  greatest  at  Point 
F.  This  is  believed  to  be  because  it  has  the  longest  period  of  sunshine. 
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Figure  6.  Computed  upward  heat  flux  for  variable  vapor  pressure, 
e,  (solid  curve;  and  for  tero  vapor  pressure  (dashed  curve). 
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The  computed  maximum  upward  heat  flux  at  all  three  sites  is 
around  0.9  cal  cm'^mln'^ , tdiich  has  the  same  order  of  magnitude  as  the 
turbulent  heat  flux  for  different  height  ranges  for  fkoe  conyeotlon 
over  flat  terrain  observed  by  Taylor  and  Swlnbank  Independently  ( see 
Priestley,  1959). 

Computations  for  temperature  and  upward  heat  flux  are  also 
made  with  two  different  assumptions  of  vapor  pressure.  The  solid  curves 
denote  the  case  for  variable  vapor  pressure  tdilch  is  based  on  observations 
made  elsewhere  for  terrain  with  vegetation.  The  dashed  cxirves  denote  the 
case  for  zero  vapor  pressure.  It  Is  seen  that  dry  atmosphere  yields  low- 
er temperatures  for  both  day  and  nighttime  because  of  decreased  long- wave 
radiation  from  the  atmosphere  and  thus  more  heat  loss  to  space.  Indeed 
this  Is  a purely  Idealistic  case.  The  solid  curve  may  be  considered  a 
case  of  high  vapor  pressure.  Other  conditioas  being  equal,  the  variation 
of  T and  heat  flux  due  to  change  of  vapor  pressure  may  be  bounded  by  these 
two  curves,  for  normal  conditions,  the  curves  should  be  closer  to  the  so- 
lid curves  for  both  temperature  and  heat  flux. 

In  view  of  the  computed  results  whown  in  figure  3.  the  eddy 
diffuslvlty  of  heat,  X , is  important  to  the  determination  of  ground 
temperature.  In  this  study  prescribed  values  are  used.  However  X should 
be  determined  ftaa  the  circulation  model  and  solved  from  the  goverming 
set  of  hydro-thezmodynamic  eqviatlons  In  the  atmospheric  boundary  layer. 

Air  temperature  is  also  a very  important  factor  in  determining  ground 
temperature.  Obtaining  air  temperature  by  solving  the  heat  equation  of 
air,  which  Includes  eddy  diffusion  and  long-wave  radiation  terms  for  air, 
would  be  very  useful. 

In  conclusion  the  computed  results  are  very  encouraging,  fhrther 


1 % 


iaproyement  can  be  made  for  application  to  mountain-valley  terrain.  The 
purpose  of  this  study  is  to  study  3-<^inensional  meso-scale  circulations 

} 

I in  typical  uneven  topography.  An  attempt  is  being  made  to  apply  this  idea 

[ to  mountain- valley  terrain.  Results  are  still  some  time  to  come. 
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